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Coherent control of quantum states is at the heart of implementing solid-state quantum pro- 
cessors and testing quantum mechanics at the macroscopic level. Despite significant progress 
made in recent years in controlling single- and bi-partite quantum systems, coherent control 
of quantum wave function in multipartite systems involving artificial solid-state qubits has 
been hampered due to the relatively short decoherence time and lacking of precise control 
methods. Here we report the creation and coherent manipulation of quantum states in a tri- 
partite quantum system, which is formed by a superconducting qubit coupled to two micro- 
scopic two-level systems (TLSs). The avoided crossings in the system's energy-level spectrum 
due to the qubit-TLS interaction act as tunable quantum beam splitters of wave functions. 
Our result shows that the Landau-Zener-Stuckelberg interference has great potential in the 
precise control of the quantum states in the tripartite system. 
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As one of three major forms of superconducting qubits ^0 a flux-biased superconducting 
phase qubit 201 consists of a superconducting loop with inductance L interrupted by a Josephson 
junction (Fig. la). The superconducting phase difference (p across the junction serves as the 
quantum variable of coordinate. When biased close to the critical current I , the qubit can be 
thought of as a tunable artificial atom with discrete energy levels that exist in a potential energy 
landscape determined by the circuit design parameters and bias (Fig. lb). The ground state |0) and 
the first excited state 1 1) are usually chosen as the computational basis states of the phase qubit. The 
energy difference between |1) and |0), u w , decreases with flux bias. A TLS is phenomenologically 
understood to be an atom or a small group of atoms tunneling between two lattice configurations 
inside the Josephson tunnel barrier, with different wave functions \L) and \R) corresponding to 
different critical current (Fig. lc). Under the interaction picture of the qubit-TLS system, the state 
of the TLS can be expressed in terms of the eigenenergy basis with \g) (the ground state) and |e) 
(the excited state). When the energy difference between |e) and \g), Hutls = E e — E g , is close 
to hui io (fa = h/2ix where h is Planck's constant.), coupling between the phase qubit and the TLS 
becomes significant, which could result in increased decoherence ^E. On the other hand, one can 
exploit strong qubit-TLS coupling for demonstrating coherent macroscopic quantum phenomena 
and/or quantum information processing^. For instance, recently a tetrapartite system formed by 
two qubits, one cavity and one TLS has been studied 5 . However, although multipartite spectral 
property and vacuum Rabi oscillation have been observed coherent manipulation of the quantum 
states of the whole system has not been demonstrated yet. 

In our experiments, we use two TLSs near 16.5 GHz to form a hybrid tripartite ^Hl phase 
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qubit-TLS system and demonstrate Landau-Zener-Stiickelberg (LZS) interference in such tripartite 
system. The avoided crossings due to the qubit-TLS interaction act as tunable quantum beam 
splitters of wave functions, with which we could precisely control the quantum states of the system. 



Results 



Experimental results of LZS interference. Fig. Id shows the measured spectroscopy of a phase 
qubit. The spectroscopy data clearly show two avoided crossings resulting from qubit-TLS cou- 
pling. Since after the application of the 7r-pulse the system has absorbed exactly one microwave 
photon and the subsequent steps of state manipulation are accomplished in the absence of the mi- 
crowave, conservation of energy guarantees that one and only one of the qubit, TLS1 and TLS2, 
can be coherently transferred to its excited state. Thus only {|l<7i<72)> |0eig 2 ), |0(?ie 2 )} as marked 
in Fig. Id are involved in the dynamics of the system. Notice that these three basis states form a 
generalized W state \ip) = a\lgig2) + /3|0eig 2 ) + 7|0<?ie 2 ), which preserves entanglement 
between the remaining bipartite system even when one of the qubit is lost and has been recognized 
as an important resource in quantum information science 13 . The system's effective Hamiltonian 
can be written as 



H = h 
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where A x (A 2 ) is the coupling strength between the qubit and TLS1 (TLS2). u T lsi (^tls2) is the 
resonant frequency of TLS1 (TLS2). Uio(t) = coioac ~ s$(t), with wio,d c being the initial energy 
detuning controlled by the dc flux bias line (i.e., the second platform holds in the dc flux bias line), 



s = \du}io(&) / d&\ being the diabatic energy-level slope of state llg^), and <&(£) being the time 
dependent flux bias (Fig. la). 

In our experiment, coherent quantum control of multiple qubits is realized with Landau- 
Zener (LZ) transition. When the system is swept through the avoided crossing, the asymptotic 
probability of transmission is exp (— 2tv^-J , where hv = dE/dt denotes the rate of the energy 
spacing change for noninteracting levels, and 2ftA is the minimum energy gap. It ranges from to 
1, depending on the ratio of A and v. The avoided crossing serves as a beam splitter that splits the 
initial state into a coherent superposition of two states These two states evolve independently 
in time while a relative phase is accumulated causing interference after sweeping back and forth 
through the avoided crossing. Such LZS interference has been observed recently in superconduct- 
ing qubits EHU However, in these experiments the avoided crossings of the single qubit energy 
spectrum are used, and microwaves, whose phase is difficult to control, are applied to drive the 
system through the avoided crossing consecutively to manipulate the qubit state. Here we use a tri- 
angular bias waveform with width shorter than the qubit's decoherence time to coherently control 
the quantum state of the tripartite system. The use of a triangular waveform, with a time resolution 
of 0.1 ns, ensures precise control of the flux bias sweep at a constant rate and thus the quantum 
state. The qubit is initially prepared in \0gig 2 ). A resonant microwave 7r-pulse is applied to coher- 
ently transfer the qubit to \ lgig 2 ). A triangular flux bias $(£) with variable width T and amplitude 

$LZS 

^x&t, t < T/2 

(2) 

2$ iZ5 -^t, T/2<t<T, 
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is then applied immediately to the phase qubit to induce LZ transitions (Fig. 2d). This is followed 
by a short readout pulse (about 5 ns) to determine the probability of finding the qubit in the state 
|1), i.e., the system in the state \lgig 2 ). 

Fig. 2a shows the measured population of |1) as a function of T and &lzs- On the top part 
of the plot, the amplitude is so small that the state could not reach the first avoided crossing Mi. 
Therefore, no LZ transition could occur and only trivial monotonic behavior is observed. When 
the amplitude is large enough to reach M 1? the emerging interference pattern can be qualitatively 
divided into three regions with remarkably different fringe patterns. 

Quantitative comparison with the model. To quantitatively model the data, we calculate the 
probability to return to the initial state Pi by considering the action of the unitary operations on 
the initially prepared state. Neglecting relaxation and dephasing, we find 

Pi = (I - Plzi) 2 + Plzi(1 ~ Plz2) 2 + Plz\Plz2 
-2P LZ1 (1 - P LZl ){l - P LZ2 ) cos(0 7 + 29 S1 - 29 S2 ) 

(3) 

-2P 2 LZ1 P LZ2 {1 - P LZ2 ) cos(0i, + 26 S2 ) 

+2P LZ1 (1 - P LZ1 )P LZ2 cos(07 + 6 n + 29 S1 ), 
where Plz% (i = 1, 2) is the Landau-Zener transition probability at the ith avoided crossing Mi, 
and Qi and Q u are the phases accumulated in region I and II, respectively (Fig. 2b). The phase 
jump 9si = Osi — tt/2 (i = 1,2) at the ith avoided crossing is due to the Stokes phase Q Si 
which depends on the adiabaticity parameter rji = Af/v in the form 6 S i = tt/4 + ^(lnr^ — 
1) + argT(l — irji), where V is the Gamma function. In the adiabatic limit 6 S — > 0, while in 
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the sudden limit 9 S = ir/4. In order to give a clear physical picture, hereafter we adopt the 
terminology of optics to discuss the phenomenon and its mechanism. First of all we define two 
characteristic sweeping rates of v x and v 2 from 27rA^/z/j = 1 (i = 1,2). From the spectroscopy 
data, we have A 1 /2n = 10 MHz and A 2 /2n = 32 MHz, thus u 1 /2n = 3.94 x 1(T 3 GHz/ns and 
u 2 /2ti = 4.04 x 10~ 2 GHz/ns, respectively. These lines of constant sweeping rate characteristic to 
the system are marked as oblique dotted lines in Fig. 2a. The avoided crossings Mi (M 2 ) can be 
viewed as wave function splitters with controllable transmission coefficients set by the sweeping 
rate v. v\ and v 2 thereby define three regions in the T — &lzs parameter plane that contain all 
main features of the measured interference patterns: 

(I) v ~ ui and v z/ 2 : M 1 acts as a beam splitter and M 2 acts as a total reflection mirror, 
i.e., Plzi ^1/2 and Plz2 — 0. In this case, equation ^ can be simplified as 

P 1 = 1 - 2P LZ1 (1 - P LZ1 )[1 + 008(0! + 26 S1 - 29 sz)}. (4) 

Apparently, only path #1 and path #2 contribute to the interference. The phase accumulated in 
region J can be expressed as 

9j= f T [uJi{t) - u 2 {t)]dt, (5) 

JO 

where Ui(t) (i = 1, 2) denotes the energy frequency corresponding to path = 1, 2). It is easy 
to find that Pi is maximized (constructive interference) in the condition 

0totai = 0j + 29 S1 - 29 S2 = (2n + l)vr, (n = 0, 1, 2, • • •), (6) 
from which we can obtain the analytical expression for the positions of constructive interference 
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fringes 

\s<$> LZS (l - T + 2~6 S i = (2n + 1)tt, S 1 < s$ LZS < 5 2 



' T + 26 S i = {2n + l)vr, 5 X < s<S> LZS < 6 2 

(7) 

^ >; . 7s + 2 (1 - J^) 5 l2 T + 2(08! - 9 S2 ) = (2n + 1)tt, s$ lzs > S 2 



where ^ = u; 1Q)(fc - c^tzsi, ^2 = ^io,dc - u TLS2 , and 5 i2 = w TL5 i - u TLS2 . 

In Fig. 2b we show the calculated constructive interference strips which agree well with the 
experimental result. Especially, in the limit of s$lzs ^> 82, $12, equation (7) can be simplified as 

<f la T + 2(0 51 -0 S2 ) = (2n+l)7r. (8) 

Intuitively, this result is straightforward to understand since in the large amplitude limit the accu- 
mulated phase 9 j is two times the area of a rectangle with length T/2 and width ojtlsi ~ ^tls2- 

(II) v ~ v 2 and u 3> V\\ Mi acts as a total transmission mirror and M 2 acts as abeam splitter, 
i.e., Plzi — 1 and Plz2 — 1/2. In this case, equation ^ can be simplified as 

Px = 1 - 2P LZ2 (1 - P LZ2 )[1 + cos(e n + 26 S2 )]. (9) 

Only path #2 and path #3 contribute to the interference. Using the same method in dealing with 
the region I, we obtain the analytical formula governing the positions of constructive interference 
fringes: 

\s§ LZS (l - ) T + 26 S2 = (2n + (10) 

* \ s<P LZS J 

As shown in Fig. 2c, the positions of the constructive interference fringes obtained from equation 



( 10) agree with experimental results very well. Similarly, in the limit s$lzs S> #2, equation ( 10) 



has the simple form 

l -s§ LZS T + 29 S2 = (2n + 1)tt, (11) 

which is also readily understood because in the large amplitude limit the accumulated phase 9 n is 
two times the area of a triangle with base-length T/2 and height s& LZ s- 

(III) V\ < v < v 2 . This region is more interesting and complex. Here, Mi acts as a beam 
splitter while M 2 can act either as a beam splitter or a total reflection mirror. This effect cannot 
be described by the asymptotic Landau-Zener formula because in this region LZS interference 
occurs only in a relatively small range around the avoided crossings. Since the analytical solution 
is extremely complicated which does not provide clear intuition about the underlying physics, 
we use a numerically calculated LZ transition probability P LZ corresponding to the transmission 
coefficient of Mi and M 2 for comparison with the experimental data. We find that for certain 
sweeping rates, LZ transition probability resulting from M 2 is quite low. Therefore, M 2 can be 
treated as a total reflection mirror while Mi is still acting as a good beam splitter. The interference 
fringes generated by M 2 thus disappear (the fringes tend to fade out) and the interference fringes 
generated by Mi dominate, displaying as a chain of 'hot spots' marked by the circles in Fig. 2a. 

When both Mi and M 2 can be treated as beam splitters, all three paths (#1, #2, and #3) 
contribute to the interference. According to equation ([3]), Pi is maximized in the condition 

6i + 2(0 S i - 9 S2 ) = (2m + 1)tt, (m = 0, 1, 2 • • •) 

(12) 

On + 29 S2 = (2n 2 + l)vr, (n 2 = 0, 1, 2 • • •) 
It is noted that under this condition the term (9i + 9n + 29 si) in equation ^ equals 2mr. Consid- 
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ering different weights in each path, it is more convenient to obtain a theoretical prediction from 
a numerical simulation. Here we utilize the Bloch equation to describe the time evolution of the 
density operator of the tripartite system: 

P = ~[H, P ]-r\p], (13) 

where T[p] includes the effects of energy relaxation. Fig. 3a shows the calculated population of |1) 
as a function of T and $ lzs- Fig. 3b is the extracted data for different T and $ LZS . The agreement 
between the theoretical and experimental results is remarkable. In order to better understand the 
origin of the 'hot spots', we also plot the probabilities of LZ transition as a function of the pulse 
width at fixed amplitude &lzs = 10 m$ (Fig. 3c). Notice that both LZ transition probabilities 
oscillate with T, which are quite different from the general asymptotic LZ transition probabilities. 
The transition probability at M 1 is always greater because A 1 is much smaller than A 2 . The three 
oblique dotted lines in Fig. 3a represent lines of constant sweeping rate. The 'hot spots' are located 
on these lines, where the transition probability of M 2 is a minimum. M 2 thereby acts as a total 
reflection mirror resulting in the 'hot spots' in transition probability. This feature further confirms 
that the avoided crossings play the role of quantum mechanical wave function splitters, analogous 
to continuously tunable beam splitters in optical experiments. The transmission coefficient of the 
wave function splitters (the avoided crossings) in our experiment can be varied in situ from zero 
(total reflection) to unity (total transmission) or any value in between by adjusting the duration and 
amplitude of the single triangular bias waveform used to sweep through the avoided crossings. 



Precise control of the quantum states in the tripartite system. We emphasize that the method 
of using LZS interference for precise quantum state manipulation described above is performed 
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within the decoherence time of the tripartite system which is about 140 ns. Through coherent LZ 
transition we can thus achieve a high degree of control over the quantum state of the qubit-TLS 
tripartite system. For example, one may take advantage of LZS to control the generalized W state, 
\ip) = a\lgig 2 ) + /3\0eig 2 ) + 7|0gie 2 ), evolving in the sub-space spanned by the three product 
states during the operation of sweeping flux bias. In order to quantify the generalized W state, we 



define w = 1 — J XXM — V V3) 2 , where a = a, f3, 7. In Fig. 3d, w is plotted as a function of T 

and &lzs- Note that with precise control of the flux bias sweep, the states with w — 1, which are 
generalized W states with equal probability in each of the three basis product states, are obtained 
demonstrating the effectiveness of this new method. It should be pointed out that when one of the 
three qubits is lost, the remaining two qubits are maximally entangled. 

In summary, our tripartite system includes a macroscopic object, which is relatively easy to 
control and readout, coupled to microscopic degrees of freedom that are less prone to environment 
induced decoherence and thus can be used as a hybrid qubit. The excellent agreement between 
our data and theory over the entire T — $lzs parameter plane indicates strongly that the states 
created are consistent with the generalized W states. The coherent generation and manipulation of 
generalized W states reported here demonstrates an effective new technique for the precise control 
of multipartite quantum states in solid-state qubits and/or hybrid qubid^l. 




Methods 



Experiment detail. Fig. la shows the principal circuitry of the measurement. The flux bias and 
microwave are fed through the on-chip thin film flux lines coupled inductively to the qubit. The 
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slowly varying flux bias is used to prepare the initial state of the qubit and to readout the qubit state 
after coherent state manipulation. In the first platform of the flux bias, the potential is tilted quite 
asymmetrically to ensure that the qubit is initialized in the left well. Then we increase the flux bias 
to the second platform until there are only a few energy levels including the computational basis 
states |0) and |1) in the left well. A microwave n -pulse is applied to rotate the qubit from |0) to 
|1). This is followed by a triangular waveform with adjustable width and amplitude applied to the 
fast flux bias line, which results in LZ transition. A short readout pulse of flux bias is then used to 
adiabatically reduce the well's depth so that the qubit will tunnel to the right well if it was in |1) 
or remain in the left well if it was in |0) . The flux bias is then lowered to the third platform, where 
the double well potential is symmetric, to freeze the final state in one of the wells. The state in 
the left or right well corresponds to clockwise or counterclockwise current in the loop, which can 
be distinguished by the dc-SQUID magnetometer inductively coupled to the qubit. By mapping 
the states |0) and |1) into the left and right wells respectively, the probability of finding the qubit 
in state |1) is obtained. We obtained Ti ~ 70 ns from energy relaxation measurement (Fig. SI a), 
Tr ~ 80 ns from Rabi oscillation (Fig. Sib), T 2 * ~ 60 ns from Ramsey interference fringe (Fig. 
Sic and Fig. Sid) and T 2 ~ 137 ns from Spin-echo (Fig. Sle) in the region free of qubit-TLS 
coupling. 

Hamiltonian in our tripartite system. For the coupled qubit-TLS system, the Hamiltonian can 
be written as 

2 2 

H q -TLSs = H q + 2J H TLSi + 22 H q _ TLSi . (14) 
i=l i=l 

In the two-level approximation the effective Hamiltonian of the qubit is H q = — |wiocrf, where 
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the flux bias ($) dependent energy level spacing of the qubit hLo w = E\ — E can be obtained 
numerically by solving the eigenvalues problem associated with the full Hamiltonian of the phase 
qubit 25 . The Hamiltonian of the ith TLS can be written as H TLSi = — ^u T LSi^ LS \ where hu T LSi 
is the energy level spacing of the ith TLS. The interaction Hamiltonian between the qubit and the 
ith TLS is H q -TLSi = ^Ajcr^ ® ^ LSl , where Aj is the coupling strength between the qubit and the 
ith TLS and o\ y (^^ s z l ) are the Pauli operators acting on the states of the qubit (the ith TLS). 
By adjusting the flux bias, the qubit and TLSs can be tuned into and out of resonance, effectively 
turning on and off the couplings. Below |0) and |1) (\gi) and |e^) ) are used to denote the ground 
state and excited state of the qubit (the ith TLS). In our experiment the initial state is prepared 
in the system's ground state \0gig 2 ). When the couplings between the qubit and TLSs are off, 
we use a 7r-pulse to pump the qubit to |1) (thus the system to \lgig 2 ). We then sweep the flux 
bias through the avoided crossing(s) to turn on the coupling(s) between the qubit and the TLS(s). 
Since after the application of the 7r-pulse the system has absorbed exactly one microwave photon 
and the subsequent steps of state manipulation are accomplished in the absence of the microwave, 
conservation of energy guarantees that one and only one of the qubit, TLSI and TLS 2, can be 
coherently transferred to its excited state. Therefore, states with only one of the three subsystems 
in excited state, \lgig 2 ), \0eig 2 ), and \0gie 2 ), are relevant in discussing the subsequent coherent 



dynamics of the system. In the subspace spanned by these three basis states the Hamiltonian ( |T4| ) 
can be written explicitly as Hamiltonian ([T]) in the main text. 

Unitary operation in our tripartite system. We use transfer matrix method l 16 l 22 l to obtain the 
probability of finding the system in \lgig 2 ) at the end of triangular pulse. We use |a)=[l, 0, 0] T , 
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|6)=[0, 1, 0] T and |c)=[0, 0, 1] T to denote the instantaneous eigenstates of the time-dependent Hamil- 



tonian (14), as shown in Fig. S2. It is noted that at the initial flux bias point, which is far from 
the avoided crossings, the system is in |a)=|lgig 2 ). At the crossing times t — ii and t = t 2 , the 
incoming and outgoing states are connected by the transfer matrix: 



cos(0i/2)exp(-i£ S i) isin(0i/2) 
isin(0i/2) cos(^i/2) exp(i0 S i) 



\ 



and 



V 

/ 



o 







(15) 







/ 



V 



(16) 



/ 



1 

cos(0 2 /2) exp(-i0<j 2 ) isin(6 2 /2) 
zsin(fl 2 /2) cos(# 2 /2) exp(i^5 2 ) 

respectively. Here sin 2 (#j/2) = PLZi(i = 1, 2) is the Landau-Zener transition probability at the ith 
avoided crossing. 6 Si = 8si — vr/2, where 8si is the Stokes phase whose value depends on the 
adiabaticity parameter rji = A?/v in the form of 6 Si = 7r/4 + r7j(ln^j — 1) +argT(l — irji), where 
T is the Gamma function. In the adiabatic limit 9 S — > 0, and in the sudden limit 8 S = tc/A. At 
crossing times t = t 3 and t = t^, we have U 3 = U 2 and U 4 = Ui, respectively. The outgoing state 
at t = U and the incoming state at t = = 0, 1, 2, 3, 4) is thus connected by the propagator 

ex P (-tjl i+1 u a (t)dt) 

exp (-i u b {t)dt) 

exp(-*/// +1 u; c (t)^ 

where Wi(t) is the energy level spacing frequency of = a, b, c) at time t. The net effect of a 



i+XA 



V 



(17) 
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triangular pulse is to cause the state vector to evolve according to the unitary transformation 

U = f>54f>4f>43f>3f>32f> 2 f>2lf>if>i0 (18) 

The probability of finding the system remaining at the initial state is P\ — \ (\g1g2\U\lg1g2) | 2 - Its 

concrete form is equation where 61 = f^[u a (t) — u b (t))dt and 6 n = J^[u b (t) — co c (t)]dt are 

the relative phases accumulated in region I and II as shown in Fig. S2, respectively. The LZS in 

our experiment can be viewed as interferences among the three paths, which are labeled with #1, 

#2, and #3, starting from the same initial state: 

path #1: \a)(t < h) -»■ \a)(ti < t < t 4 ) -»■ \a)(U < t < T) 

path #2: \a)(t < h) -»■ |6)(ti < t < t A ) ->• \a)(t 4 < t < T) 

path #3: \a)(t < h) ->■ |6)(*i < t < i 2 ) -» |c> (t 2 < t < t 3 ) -> |6)(t 3 < t < t 4 ) ->■ |a)(*4 < * < 
T) 

Denoting Wj(t) as the energy level spacing frequency corresponding to path i(i = 1, 2, 3), then 9i 
and 9n have the forms 8j = / T [u;i(t) — u>2(t)]dt and On = Solvit) — u)^{t)\dt, respectively. 

Numerical simulation of LZS interference in bipartite qubit-TLS system. For the bipartite 
qubit-TLS system discussed here, the qubit is coupled only to a single TLS. The quantum dynamics 
of the system, including the effects of dissipation, is described by the Bloch equation of the time 
evolution of the density operator: 



p=~[H b ,p]-T\p], (19) 



where H b = h 



f wio(t) A ^ 



with u;io(t) = tdio,dc — vt, v = 2s$lzs/T is the energy sweep- 

y A ootls 1 

ing rate and A is the qubit-TLS coupling strength. The second term T[p] describes the relaxation 
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process to the ground state \0g) and dephasing process phenomenologically. In a concrete expres- 



sion, equation (19) can be written as: (for ease of discussion, we relabel \ lg) and |0e) as \a) and 



| b), respectively) 



p aa = -iA(p ba - p ab ) - T a p aa , 

Pbb = i^iPba ~ Pab) - ^bPbb, 

p ab = ~iA(p bb - Paa) - i(Uio(t) - U T Ls)Pab 



(20) 



^YabPab) 

with p ba = p* ab - Here T a (a = a, b) is the relaxation rate from state |a) to the ground state \0g). 
The decoherence rate ^ ab = (T a + T b )/2 + ^ de P h ^ includes contributions from both relaxation and 
dephasing. Fig. S3a and Fig. S3b give the numerically simulated LZS interference pattern for 
the qubit coupled with the first TLS and second TLS, respectively. To calculate the transmission 
coefficient of M f (i = 1,2), i.e., the Landau-Zener tunneling probability P LZ , as shown in Fig. 3c 
in the main text, we cannot directly use the asymptotic Landau-Zener formula which is based on 
sweeping the system across the avoided crossing from negative to positive infinities. In contrast, 
in our experiment the LZS occurs near the avoided crossings. Therefore, our numerical results are 
obtained by solving the Bloch equations directly. 



Numerical simulation of LZS interference in tripartite qubit-TLS system. For the tripartite 
qubit-TLS system discussed below, the qubit is coupled resonantly to two TLSs (TLS1 and TLS2) 
with different excited state energies Hcutlsi an d Tiujtls2- The Hamiltonian in the basis of \lgig 2 ), 
|0ei<72)> |0gie 2 ) is Hamiltonian ([TJ) in the main text. The Bloch equations that govern the evolution 
of the density operator can be written as: (For simplicity, we relabel \lgig 2 ), \0eig 2 ), |0gie 2 ) as 
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\a), \b), | c), respectively.) 

Paa = ~iAl(Pba ~ Pab) ~ ^(Pca ~ Pac) ~ ^aPaa 
p bb = tA^Pba ~ Pab) ~ FbPbb 

Pec = ?A 2 (p ca - p ac ) - Y c p cc 

Pab = -i(u) W (t) - U T LSl)Pab ~ ^l{Pbb ~ Paa) 

- iA 2 pcb - labpab W 
Pac = -i(Ul (t) - UJTLS2)Pac ~ «A 2 (p cc ~ Paa) 

- tAxPbc - lacPac 

p hc = -i(u TLS i - u TLS2 )pbc - ^Aip ac 

+ iA 2 p ha - IbcPbc 

V 

where the diagonal elements pa are the populations, off-diagonal elements p^i ^ j) describe 
coherence, and 7^ = (r, + Tj)/2 + ^ de v h "> are the rates of decoherence. The remaining three 
elements' equations are determined by p*j = pji. The numerically simulated LZS interference 
pattern is shown in Fig. 3a in the main text, which agrees with experimental result excellently. 
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Figure 1 Qubit circuit and experimental procedure, a, Schematic of the qubit cir- 
cuitry. Josephson junctions AI/AIOx/AI are denoted by the X symbols. The flux bias, 
microwave and readout dc-SQUID are inductively coupled to the qubit with inductance 
L w 770 pH, capacitance C « 240 fF and critical current I w 1.4 fxA. b, Principle of the 
operation and measurement of the phase qubit. The two lowest eigenstates |0) and |l) 
form the qubit with transition frequency u 10 which can be adjusted by changing the flux 
bias. A microwave pulse is used to manipulate the qubit state and readout pulse then 
lower the potential energy barrier to perform a fast single-shot readout, c, Schematic of a 
two-level state located inside the insulating tunnel barrier of a Josephson junction and its 
eigenstates in different bases, d, Spectroscopy of the coupled qubit-TLS system with cor- 
responding quantum states labeled. Two avoided crossings centered at u TLS i and u TLS 2 
are observed. 

Figure 2 LZS interference in a phase qubit coupled to two TLSs. a, The population 
of |l) measured immediately (a few ns) after the triangular flux pulse is plotted as a func- 
tion of the width and amplitude of the triangular flux bias waveform. The oblique dotted 
lines are lines of constant characteristic sweeping rate v x and v 2 defined in the text. The 
white circles mark the 'hot spots', where the interference fringes generated by M 2 tend 
to fade out and the interference fringes generated by Mi dominate, b and c, Analytically 
calculated constructive interference strips in region I and II, respectively. The horizontal 
and vertical dotted lines indicate the corresponding locations of interference strips, b and 
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c have the same axis labels as a. d, Schematic of generating LZS interference with tun- 
able beam splitters in a phase qubit coupled to two TLSs. Mi and M 2 correspond to the 
TLSs with smaller and larger avoided crossings in Fig. "Id, respectively. 



Figure 3 Numerically simulated LZS interference pattern and control of a general- 
ized W state in a phase qubit coupled to two TLSs. a, The numerically simulated popu- 
lation of |l) after the triangular flux pulse is plotted as a function of the width and amplitude 
of the triangular flux bias. The horizontal dotted line indicates the location of $ L zs = 10 
m$o and the vertical dotted lines indicate the locations of 'hot spots' at § LZS = 10 m$ . 
The oblique dotted lines are lines of constant sweeping rate. The parameters used are 
determined experimentally: u 014c /2n = 16.747 GHz, \s\ = ||f | = 0.0404 GHz/m$ , 
wtl5i/2tt = 16.590 GHz, co TLS2 /2n = 16.510 GHz, A 1 /2ir = 10 MHz, A 2 /2tt = 32 MHz, 
T lgig2 = (70 ns)- 1 , r 0ei92 = r 09ie2 = (146 ns)' 1 , 7 (***) = (45 ns)- 1 . b, The upper panel 
shows the dependence of population of |l) on $ LZ5 at T = 20 ns, 40 ns, 60 ns, respec- 
tively. The lower panel shows the dependence of population of |l) on T at § LZ s = 3.6 
m$o, 7.2 m$o, 10.8 m$ , respectively. The circles represent the experimental data and the 
lines from the theory, c, LZ transition probabilities of M 1 (blue line) and M 2 (red line) at 
§lzs = 10 m$ as a function of pulse width. They are quite different from the asymptotic 
LZ transition probabilities (blue dotted line and red dotted line), d, The resulting w as a 
function of T and § L zs- 
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